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ABSTRACT 

In  this  thesis,  an  effect  of  turbulent  gas  flow  on  the  glow  dis- 
charge-to-arc transition  is  investigated.  This  is  accomplished  by 
solving  the  continuity  equation  for  charged  particles  by  two  different 
numerical  processes.  They  are  a  linearized  solution  and  a  non- linear 
numerical  solution.  The  numerical  results  using  an  axi-symmetric 
cylindrical  geometry  show  the  charge  particle  density  profile  with 
respect  to  radial  distance,  time  and  diffusion  coefficient. 

In  nitrogen  at  a  time  equal  to  4x  ,  the  radius  of  the  streamer  in 

near  turbulent  conditions  (p  =  50)  was  found  to  ba  2.5  times  larger  than 

4 
in  laminar  conditions  (p  =  10  ).   i  is  the  characteristic  recombination 

time  and  p  is  the  ratio  of  the  characteristic  diffusion  time  divided  by 

V 

The  effective  increase  in  the  diffusion  coefficient  between  laminar 
and  turbulent  gas  flow  causes  the  charge  density  profile  to  expand  more 
radially  outward.  This  reduces  the  charge  density  within  the  original 
streamer  volume  tube  causing  a  reduction  in  the  conductivity  of  the 
streamer.  This  reduction  in  conductivity  may  possibly  lead  to  a  delay 
in  the  development  of  a  breakdown  streamer.  As  a  result,  the  gas  may  be 
raised  to  higher  current  levels  before  breakdown  occurs. 
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I.   INTRODUCTION 

The  belief  that  a  turbulent  flow  of  an  electrically  excited  lasing 
medium  leads  to  an  increase  in  the  glow-discharge  stability  has  been 
demonstrated  for  a  number  of  years  [Ref.  1].  However,  few  experiments 
that  yield  definitive  results  have  been  reported.  The  increase  in  glow 
discharge  stability  has  been  attributed  to  the  dramatic  increase  in  the 
diffusion  transport  mechanism  and  thermal  conductivity  of  the  gas 
between  laminar  and  turbulent  flow.  This  may  rapidly  dilute  the  local 
charge  carrier  density  in  the  gas  and  thus  delay  the  onset  of  the 
streamer-to-arc  transition.  Previous  experimental  and  theoretical 
analyses  have  centered  around  the  prevention  of  the  breakdown  discharge; 
that  is  the  glow-discharge-to-arc  transition  [Ref.  2].  This  is  done  by 
observing  the  current-to-voltage  relationship  of  a  gas  with  and  without 
turbulent  flow.  These  investigations  generally  do  not  attempt  to  answer 
the  question  of  what  is  the  effect  of  turbulent  flow  on  the  onset  of  the 
arc.  This  work  proposes  to  analyze  the  events  that  reduce  the  charge 
carrier  density  in  the  arc  channel.  With  this  in  mind,  a  model  is 
presented  which  describes  the  onset  and  delay  mechanisms  involved  in  the 
turbulent  flow,  glow-discharge-to-arc  transition. 

Even  though  we  pretend  to  study  and  model  the  effects  of  turbulent 
flow  on  the  glow-discharge-to-arc  transition  in  an  electrically  pumped 
laser,  the  information  provided  should  be  applicable  to  the  onset  of 
many  spark  or  arc  situations. 
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The  narrative  is  written  at  a  level  where  a  graduate  student  with 
fundamental  skills  in  mechanics  and  an  introduction  to  basic  atomic 
processes  will  be  able  to  understand  the  subtleties  of  discharges.  With 
such  a  limited  background  in  mind,  a  physical  model  of  the  ionization 
and  sparking  process  is  described.  The  process  of  ionization  and  de- 
ionization  are  also  reviewed. 

The  presentation  will  begin  with  the  modern  understanding  of  ioni- 
zation in  a  gas.  Then,  theories  which  pertain  to  the  development  of  the 
transition  to  an  arc  will  be  discussed  and  experimental  results  per- 
taining to  the  turbulent  flow  arc  transition  will  be  introduced.  After- 
wards, an  analysis  of  the  streamer  decay  will  be  developed.  Finally, 
the  results  of  a  numerical  analysis  will  be  used  to  substantiate  the 
delay  mechanism  theorized  with  the  glow-discharge-to-arc  transition  in  a 
turbulent  flow  situation. 

The  following  is  a  description  of  an  arc  in  a  gaseous  medium.  It 
relates  the  voltage-current  relationship  in  a  gas  and  tells  in  a  quali- 
tative manner  how  an  arc  occurs  and  its  relation  to  the  glow-discharge. 

"In  the  typ-'cal  voltage-current  characteristic  of  a  gas  the  last 
region  of  discharge,  that  is,  the  one  having  the  highest  current,  is 
called  the  'arc1.  In  this  region,  the  current  is  about  two  orders  of 
magnitude  higher  than  in  the  glow,  whereas  the  voltage  is  roughly  an 
order  of  magnitude  lower. 

The  arc  is  probably  the  form  of  the  discharge  least  understood 
theoretically,  although  it  has  been  used  for  a  long  time  in  a  number 
of  applications. 

The  glow-to-arc  transition  occurs  when  the  abnormal  glow  has 
acquired  a  critical  value  of  the  voltage  and  when  the  power  supply  has 
an  adequately  low  internal  resistance  capable  of  sustaining  the  high 
currents  required  for  the  arc. 
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Arcs  have  been  known  since  the  early  days  of  electricity  because 
of  their  immense  light  intensity,  which  has  been  and  still  is  used  for 
lighting  purposes.  Because  of  the  arc's  high  temperature,  extensive 
research  has  been  recently  conducted  to  produce  the  temperatures 
needed  to  trigger  thermonuclear  fusion. 

In  the  arc  the  current  density  is  much  higher  than  the  glow.  The 
increased  light  intensity  makes  it  very  difficult  to  recognize  the 
various  regions  within  the  arc.  At  atmospheric  pressure  in  air,  the 
arc  seems  to  have  a  very  bright  core  with  less  bright  surroundings 
which  resemble  a  hot  flame.  The  core  itself  has  no  discernable  struc- 
ture. An  arc  keeps  changing  position  because  of  the  continuous 
heating  that  leads  to  convection  of  the  hot  gas  upward.  A  horizontal 
arc  will  therefore  have  an  arched  shape,  from  which  it  derived  it's 
name. 

The  arc  is  a  very  hot  plasma.  Usually,  arcs  are  produced  in  the 
near  atmospheric  range  and  therefore  the,££harged-particL2  density  in 
their  cores  will  range  from  10    to  10    electrons/cm  ."  [Ref.  3]. 

Many  familiar  discharge  events  surround  us  every  day.  The  light 
emitted  from  a  fluorescent  lamp  is  a  glow-discharge  type  event  which 
occurs  at  a  reduced  pressure.  Lightning  is  a  large  scale  ele:tric 
discharge  at  atmospheric  pressure.  Glow  discharges  are  often  usud  to 
excite  electrons  to  high  energy  states  in  electrically  pumped  qas  lasers. 
A  breakdown  or  arc  within  the  laser  can  cause  structural  damage.  More 
importantly,  arcing  destroys  the  population  inversion  needed  to  acnieve 
lasing  action.  To  obtain  the  maximum  power  output  from  this  type  of 
laser  requires  operating  the  glow  discharge  'pump'  very  near  arcing 
conditions.  With  the  addition  of  turbulent  flow  of  the  gas,  tiie  transi- 
tion from  glow  discharge  to  arc  is  delayed.  This  permits  either  a 
greater  margin  of  safety  in  laser  operation  or  higher  power  levels 
produced  by  the  laser. 
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II.  METHODS  OF  IONIZATION 

"The  process  of  liberating  an  electron  from  a  gas  with  either  the 
simultaneous  production  of  a  positive  ion  or  with  the  increase  of 
positive  ion  charge  is  called  ionization"  [Ref.  4]. 

The  two  most  important  ionization  processes  in  a  gas  are  by  electron 
collision  and  by  the  absorption  of  radiation  (photoionization).  If 
ionization  if  a  neutral  atom  is  to  uccur  by  collision  directly,  the 
impinging  electron  must  have  a  kinetic  equal  to  or  greater  than  the 
ionization  energy  of  the  particular  atom.  Atoms  may  also  be  ionized  in 
steps  by  electron  collisions  with  energy  less  than  that  needed  for 
direct  ionization.  This  is  a  less  efficient  method  of  ionization, 
however. 

H.  Hertz  found  that  a  spark  between  two  electrodes  at  a  lower  volt- 
age could  be  created  when  the  light  emission  from  another  spark  were 
permitted  to  impinge  on  it.  He  showed  that  the  ultraviolet  light  from 
the  first  spark  striking  the  cathode  of  the  second  caused  the  breakdown. 
The  ionization  energy  of  hydrogen,  for  example,  is  about  14  electron 

volts.   For  a  photon  to  ionize  hydrogen  directly  requires  that  it  have  a 

15 
frequency  of  about  3.4  •  10   Hz,  which  is  in  the  near  ultraviolet. 

This  assumes  one-hundred  percent  efficiency  and  direct  liberation  of  an 

electron  due  to  the  photon  absorption  by  the  hydrogen  atom. 

When  electrodes  are  put  into  a  gas  such  as  air  or  nitrogen  and  a 

weak  voltage  is  applied  to  their  terminals,  a  very  weak  current,  much 

less  than  a  picoampere  by  several  orders  of  magnitude  will  appear.  This 


current  will  increase  with  voltage  until  all  of  the  charge  carriers 
naturally  present  are  being  used  and  current  saturation  occurs.  This 
means  that  all  of  the  charge  carriers  move  to  the  electrodes  with  little 
loss  due  to  recombination  and  diffusion  to  the  walls.  These  charge 
carriers  are  created  naturally  by  external  radiation  such  as  cosmic 
rays.  Further  increases  in  the  applied  voltage  beyond  the  plateau 
discussed  above  will  yield  a  rapidly  increasing  current  flow  on  the 
order  of  a  picoampere.  This  zone  is  called  the  Townsend  discharge 
region  and  the  new  charge  carriers  are  being  created  by  ionizing  colli- 
sions between  fast  moving  electrons  and  neutral  atoms.  These  fast 
moving  electrons  are  accelerated  by  the  electric  field  and  upon  collid- 
ing with  a  neutral  gas  particles  may  either  excite  the  gas  or  ionize  it, 
further  increasing  the  density  of  charged  particles.  This  ionization 
will  take  place  when  the  kinetic  energy  gained  by  the  electron  between 
collisions  is  greater  than  that  required  for  ionization.  Clearly,  the 
lower  the  gas  density,  the  farther  the  electron  will  travel  between 
collisions  (lower  collision  frequency)  aid  the  larger  will  be  it's 
average  kinetic  energy. 

In  this  scheme,  the  total  number  of  available  electrons  will  in- 
crease if  the  initial  electrons  have  the  energy  required  to  ionize  the 
gas  and  a  self-sustaining  current  will  develop.  The  new  electrons  will 
gain  energy  in  the  same  manner  and  the  number  of  electrons  will  increase 
between  the  cathode  and  the  anode.  This  process  is  called  an  electron 
avalanche  due  to  the  rapid  increase  in  electron  production. 
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After  this  interval,  the  current  increases  by  orders  of  magnitude 
with  a  slight  increase  of  applied  voltage.  Here,  the  current  is  in- 
dependent of  external  ionizing  sources  and  the  current  becomes  self- 
sustaining.  Further  increases  in  the  current  flow  by  reducing  the 
external  limiting  resistor  produces  a  glow  in  the  gas.  This  is  the 
primary  region  of  interest  in  an  electrically  pumped  laser. 

Glow-discharge  is  one  of  the  oldest  known  discharge  phenomena.  It 
develops  very  easily  when  the  gas  pressure  is  less  than  one-hundred 
Torr.  Glow-discharge  will  also  appear  under  special  conditions  at  higher 
pressures.  The  diffuse  glow-discharge  region  is  well  suited  as  a  laser 
pumping  scheme.  This  is  primarily  due  to  the  relatively  good  coupling 
of  electron  energy  to  excited  atomic  states.  These  states  can  be  either 
vibrational  or  electronic  in  nature. 

When  electrons  leave  the  cathode,  they  are  accelerated  by  the 
applied  electric  field.  As  the  electrons  are  accelerated,  they  quickly 
become  too  energetic  to  recombine  with  passing  ions.  As  the  electrons' 
energy  increase,  they  ionize  neutral  atoms  by  collision.  This  process 
retards  the  motion  of  the  electrons  enough  that  they  are  now  able  to 
recombine  with  passing  ions.  The  result  of  this  two-body  recombination 
is  the  emission  of  a  photon  which  carries  off  excess  energy  in  the  form 
of  light.  In  the  glow-discharge  this  region  is  known  as  the  negative 
glow. 

Finally,  the  transition  to  an  arc  discharge  occurs  when  the  current 
is  allowed  to  increase  further.  A  highly  ionized  channel  forms  which 
collapses  the  glow.  The  current  through  this  channel  is  generally 
limited  only  by  the  external  circuit  network.  The  arc  breakdown  of  the 
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gas  is  an  irreversible  transition.  If  this  plasma  channel  is  short 
lived  and  at  atmospheric  pressures,  it  is  known  as  a  spark.  If  the 
current  is  steady,  the  degree  of  ionization  is  determined  by  the  magni- 
tude of  the  current  and  a  plasma  state  known  as  an  arc  is  developed. 

There  are  only  two  generally  accepted  breakdown  theories.  The  first 
was  developed  by  Townsend  in  the  early  1900s.  In  this  method,  succes- 
sive electron  avalanches  between  the  cathode  and  anode  develop  in  which 
every  avalanche  produces  at  least  one  following  avalanche.  A  highly 
ionized  channel  is  produced  in  which  the  conductivity  becomes  so  high 
that  the  current  through  it  is  limited  by  the  external  circuit  elements. 
This  theory  remained  unchallenged  until  the  development  of  high  speed 
oscilloscopes.  With  these  improved  observation  tools,  it  was  found  that 
the  development  times  of  the  ionized  channel  was  very  much  faster  than 
the  Townsend  theory  could  predict.  The  period  of  time  required  for  the 
irreversible  transition  is  known  as  the  formative  lag  of  breakdown. 
Formative  time  lags  on  the  order  of  ten  nanoseconds  were  first  recorded 
in  the  1920s.  These  short  times  almost  totally  precluded  the  use  of 
positive  ion  schemes  in  breakdowns  due  to  their  relatively  large  mass 
and  limited  the  role  of  the  electron  avalanches. 

The  ionization  channels  were  seen  to  be  heavily  branched  and  fila- 
mentary in  nature.  They  could  be  started  at  low  voltages  at  asymmetric 
nodes  and  could  reach  large  distances.  This  channel  is  also  known  as  a 
streamer.  If  breakdown  occurs  by  a  transition  from  a  streamer  to  a 
spark,  the  mechanism  involved  is  often  called  a  streamer-breakdown.  The 
streamers  develop  in  time  intervals  on  the  order  of  ten  nanoseconds 
[Ref.  5]. 
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The  second  method  which  causes  breakdown  to  occur  is  based  on  pho- 
tons which  ionize  the  gas.  During  the  buildup  of  the  primary  avalanche, 
excitation  of  the  gas  as  well  as  ionization  takes  place.  These  excited 
states  have  lifetimes  on  the  order  of  picoseconds.  Thus,  before  the 
avalanche  can  reach  its  full  size,  photons  will  be  emitted  from  these 
excited  states  as  they  return  to  the  ground  state. 

The  photons  will  be  emitted  isotropically  and  will  be  absorbed  at 
various  distances  from  their  origin,  depending  on  the  absorption  co- 
efficient of  the  gas.  Numerous  processes  then  take  place  which  lead  to 
photoionization.  Newly  liberated  electrons  are  then  available  at 
various  points  in  the  gas  due  to  the  initial  avalanche  tnat  :s  still 
advancing  toward  the  anode.  The  photoelectrons  will  be  accelerated  by 
the  electric  field  and  can  cause  avalanches  of  their  own.  Tl"ese  new 
avalanches  can  repeat  the  whole  process  and  thus  account  for  tie  short 
formative  times  involved  in  the  streamer  breakdown  process  [Ref.  6]. 

The  positive  ions,  due  in  part  to  their  much  larger  mass,  will  be 
virtually  left  behind  as  the  electron  avalanche  advances.  They  consti- 
tute the  highly  ionized  channel  that  is  characteristic  of  the  spark. 
Those  ions  act  to  carry  further  electrons  which  sustain  the  discharge. 

This  streamer  mechanism  of  breakdown  is  much  faster  than  that  of  the 
Townsend  mechanism.  Single  avalanches  can  produce  streamers  which 
proceed  at  a  much  shorter  time  than  the  electron  crossing  time.  This 
accounts  for  the  short  formative  time  lags  [Ref.  7]. 

Streamers  are  the  forerunners  of  the  breakdown  process.  In  a  uni- 
form electric  field,  streamer  formation  always  leads  to  breakdown.   In 
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fact,  breakdown  can  also  occur  without  the  development  of  streamers.  It 
has  been  shown  [Ref.  8]  that  there  is  a  minimum  avalanche  size,  below 
which  streamers  do  not  form. 

As  already  stated,  there  are  only  two  known  methods  which  are  re- 
sponsible for  breakdown.  Experiments  show  [Ref.  9].  that  there  is  no 
abrupt  change  that  would  signal  a  transition  from  one  type  to  the  other. 
Often  the  discharge  follows  the  Townsend  pattern  followed  by  a  streamer 
occuring  near  the  anode.  This  streamer  then  propagates  to  the  cathode 
causing  breakdown.  This  pattern  is  best  explained  by  considering  the 
cloud  of  positive  ions  which  is  left  behind  the  advancing  avalanche. 
This  distorts  the  local  electric  field  in  such  a  manner  as  to  enhance 
the  avalanche  effect. 

The  basic  difference  between  breakdown  in  a  nonuniform  field  in 
contrast  with  a  uniform  field  is  the  fact  that  streamers  which  develop 
in  a  uniform  field  always  lead  to  breakdown.  In  the  nonuniform  case, 
however,  a  variety  of  events  take  place  before  the  breakdown  occurs. 
The  first  experimenters  were  interested  in  the  glow  that  would  appear 
before  the  breakdown.  These  glows,  or  better  yet  glow-discharges,  were 
either  self-sustaining  or  transient  depending  primarily  on  the  degree  of 
electric  field  nonuniformity.  These  coronas  have  been  observed  natur- 
ally during  electrical  storms  and  as  "St.  Elmo's  Fire".  They  generate  a 
considerable  amount  of  electrical  noise  as  well. 

The  transition  to  glow  from  intermittent  streamers  is  thought  to  be 
due  to  the  development  of  a  space  charge  around  the  anode.  This  space 
charge  creates  somewhat  of  a  uniform  electric  field  in  this  region. 
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This  can  lead  to  a  breakdown  within  this  gap  which  is  responsible  for 
the  glow.  This  type  of  discharge  is  relatively  stable  due  to  a  variety 
of  mechanisms  [Ref.  10].  Such  glow  will  only  develop  in  the  presence  of 
negative  ion  formation.  Corona  glow  raises  the  arcing  threshold  voltage 
which  would  otherwise  lead  to  a  streamer-to-arc  breakdown  transition. 

When  voltage  is  applied  to  electrodes  so  as  to  create  a  highly 
nonuniform  electric  field,  the  first  observable  phenomenon  is  that  of  a 
branched  filament  already  identified  as  a  streamer.  The  more  uniform 
the  electric  field,  the  less  branching  of  the  anode  streamer.  Experi- 
ments shown  that  the  branches  never  cross  and  tend  to  repel  each  other 
[Ref.  11].  The  number  of  branches  increases  from  the  cathode  to  the 
anode.  This  increase  is  exponential.  The  anode  streamer  channels  can 
thicken  due  to  new  growth.  The  glow  from  these  streamers  comes  only 
from  the  tips  which  are  advancing.  This  is  where  the  new  ionization  is 
taking  place. 

Electron  emission  from  the  cathode  can  take  place  when  the  striking 
anode  streamers  are  of  sufficient  intensity.  The  emission  is  due  to  the 
large  electric  field  of  the  space  charge  in  front  of  the  anode  streamer 
tips.  Electrons  accelerated  toward  the  streamer  can  create  avalanches 
of  new  electrons  by  collisions  with  the  neutral  gas.  Thus  the  conduc- 
tivity is  rapidly  enhanced  and  a  jump  in  the  current  between  the  elec- 
trodes is  observed. 

A  fully  ionized  plasma  behaves  differently  than  a  partially  ionized 
in  many  ways.  One  principle  difference  is  in  the  conductivity.  It  can 
be  shown  [Ref.  12]  that  the  conductivity  of  a  fully  ionized  plasma  is 
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proportional  to  the  electron  energy  to  the  3/2  power.  However,  in  a 
partially  ionized  plasma,  the  conductivity  is  dependent  on  the  charged 
particle  density  directly. 

The  conductivity  of  a  streamer  is  not  very  high.  Therefore,  break- 
down may  not  occur  when  the  streamer  channel  reaches  the  cathode.  It  is 
the  electrons  emitted  from  the  cathode,  propagating  through  the  channel, 
that  increases  the  overall  conductivity. 

When  the  voltage  is  slowly  applied  to  the  electrodes,  a  few  streamer 
bursts  occur.  These  are  principally  due  to  random  fluctuations  in  the 
gas.  As  the  voltage  is  increased  further,  streamers  develop  more  often 
and  eventually  lose  their  fluctuating  nature.  When  the  transients  die 
away,  a  steady  glow  appears  near  the  anode  surface.  The  current  between 
the  electrodes  is  continuous  with  small  fluctuations  superimposed  on  it. 
With  increasing  voltage  and  hence  increasing  current,  the  glow  increases 
in  both  intensity  and  size  [Ref.  13].  Just  prior  to  the  onset  of  break- 
down, a  streamer  type  transient  discharge  appears  with  the  associated 
increase  in  anode  glow  brightness  and  electrical  noise.  This  is  called 
a  breakdown  streamer  whereas  the  former  is  known  as  an  onset  streamer. 

Figure  1  [Ref.  14]  shows  the  different  current-voltage  regions  which 
occur  in  dry  air  at  atmospheric  pressure.  The  electric  field  is  highly 
non-uniform  due  to  special  electrode  arrangements.  The  first  region  is 
the  current  produced  by  natural  ionization  of  the  gas  due  primarily  to 
cosmic  rays.  The  next  region  is  that  caused  by  ionization  due  to 
collisions.  This  is  followed  by  increased  current  from  the  onset 
streamers.   When  the  sharp  current  increase  slacks  off,  a  steady  glow 
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Figure  1.   Dry  Air  Voltage-Current  Diagram 
[Ref.  14] 
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forms  which  transitions  into  breakdown  streamers  and  finally  a  spark  or 
breakdown  with  increased  voltage. 

Breakdown  streamers  have  been  seen  to  have  a  slow  rise  and  longer 
decay  time  associated  with  them  when  compared  to  onset  streamers.  These 
pulses  occur  with  a  frequency  from  zero  to  three  to  four  thousand  pulses 
per  second.  Beyond  this  the  frequency  decays  but  the  streamer  duration 
increases. 

Near  a  cathode  that  has  a  highly  nonuniform  electric  field,  similar 
events  occur  which  lead  to  cathode  glow  [Ref.  15]. 
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III.   AMBIPOLAR  DIFFUSION 

Due  to  a  density  gradient  within  the  plasma  channel,  there  will  be 
an  exodous  of  mass  due  to  diffusion.  If  the  electron  density  is  low,  as 
in  a  partially  ionized  gas,  each  species  of  particles  may  be  considered 
independently.  However,  this  may  not  be  the  case  in  the  highly  ionized 
spark  channel . 

Diffusion  prevents  a  well  defined  boundary  between  the  highly  ion- 
ized channel  and  the  slightly  ionized  gas  which  surrounds  it  from 
existing  for  all  but  very  short  times.  If  the  distribution  of  electrons 
and  ions  is  not  well  behaved,  (as  in  the  presence  of  magnetic  fields), 
then  the  electrostatic  repulsion  among  like  particles  will  assist  in  the 
diffusion  process.  This  is  precisely  the  case  that  occurs  when  there  is 
a  build-up  of  either  electronic  or  ionic  space  charges.  The  accurate 
treatment  of  electrostatic  repulsion  is  beyond  the  scope  of  this  work 
and  will  not  be  treated  further. 

To  create  a  steady-state  discharge  situation,  the  need  exists  to  not 
only  ionize  the  gas,  but  also  to  de-ionize  the  gas.  The  rate  of  ioniza- 
tion must  equal  the  rate  of  de-ionization.  The  loss  of  electrons  from 
the  spark  channel  occurs  principally  in  two  ways:  diffusion  out  of  the 
system  or  recombination  of  an  ion  with  an  electron.  These  lead  to  the 
direct  de-ionization  of  the  plasma. 

Electron  attachment  to  a  neutral  atom  may  also  occur,  especially  in 
a  medium  with  electro-negative  atoms  such  as  oxygen.  The  negative 
charge  is  not  lost  to  the  system;  however,  the  large  increase  in  the 
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mass  of  the  body  reduces  the  charge's  mobility  by  several  orders  of 
magnitude.  Positive  ions  may  perform  the  same  process  which  is  called 
cluster  formation.  Since  the  mobility  of  the  ion  is  so  much  smaller 
than  the  electron  to  begin  with,  the  overall  mobility  of  the  cluster 
formation  is  not  reduced  to  the  degree  that  the  attachment  process  slows 
the  electron. 

Collisions  between  like  particles  give  rise  to  no  diffusion.  This 
situation  is  to  be  contrasted  with  the  case  of  ions  colliding  with 
neutral  atoms.  The  final  velocity  of  the  neutral  is  of  no  concern  and 
the  ion  random-walks  away  (diffuses)  from  its  initial  position.  Ion-ion 
collisions  have  a  balance  in  each  collision.  That  is,  for  each  ion  that 
moves  out,  the  other  moves  in  as  a  result  of  the  collision. 

When  an  electron  and  an  ion  collide,  however,  the  situation  is 
different.  Unlike  particle  collisions  give  rise  to  diffusion.  The 
physical  picture  is  difficult  due  to  the  disparity  in  mass  between  the 
electron  and  ion.  The  ions  remain  relatively  stationary  and  the  elec- 
trons do  the  random  walk  [Ref.  16]. 

Particles  that  are  able  to  diffuse  to  the  wall  of  the  container  are 
able  to  recombine  and  are  thus  lost  to  the  system.  Therefore,  within 
the  gas  the  electrons  can  thermally  diffuse  to  the  walls,  recombine 
within  the  plasma,  create  negative  ions,  or  undergo  ambipolar  diffusion. 

Under  the  influence  of  thermal  forces,  the  electrons'  smaller  mass 
than  the  ions'  means  thay  can  diffuse  faster  than  the  ions.  In  the 
process,  the  electrons  leave  behind  an  excess  number  of  ions  which 
creates  a  local  space  charge  and  thus  an  electric  field.  This  micro- 
scopic field  will  then  tend  to  retard  the  motion  of  further  electrons 
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and  propel  the  ions.  The  two  types  of  charge  carriers  will  diffuse  as 
one.  This  is  called  ambipolar  diffusion  [Ref.  17].  The  effect  of  the 
ambipolar  electric  field  is  to  enhance  the  diffusion  of  ions  by  a  factor 
of  two,  but  the  diffusion  rate  of  the  electrons  and  ions  together  is 
primarily  controlled  by  the  slower  species. 

Electrons  which  cannot  propagate  to  the  electrodes  will  be  lost 
through  recombination,  diffusion  or  attachment.  The  charge  carriers 
will  diffuse  throughout  the  medium  if  there  exists  a  density  gradient  in 
much  the  same  manner  as  neutral  atoms  if  there  are  no  forces  to  confine 
them.  When  they  reach  the  wall,  they  recombine  and  are  lost  to  the  gas. 

The  current  density  due  to  diffusion  is  [Ref.  18]. 
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The  rate  of  change  of  particle  density  due  to  diffusion  is 
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If  we  assume  that  all  ionization  mechanisms  have  stopped,  then  the  rate 
of  change  of  electron  density  due  to  diffusion  is: 
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[Ref.  19]     (3) 


In  many  gaseous  discharges  the  charged  particles  are  contained 
within  a  cylindrical  channel.  This  is  the  exemplified  in  the  fluor- 
escent light  tibes.  If  the  channel  is  much  longer  than  it  is  wide,  then 
the  end  effects  may  be  neglected.  In  addition,  the  radius  of  the 
cylinder  must  be  larger  than  the  Debye  length.  The  analysis  of  the 
effects  of  diffusion  and  recombination  are  therefore  done  in  a  cylin- 
drical geometry.  Moreover,  the  diffusion  coefficient  (D)  is  simply 
treated  as  a  constant  in  equation  (3).  Note  that  the  net  current  in 
ambipolar  diffjsion  is  zero  as  the  electron  and  ion  diffuse  together. 
However,  the  ambipolar  particle  density  change  is  driven  by  density 
gradients  rather  than  electro-static  forces. 
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IV.   RECOMBINATION 

The  term  recombination  refers  to  a  charge  neutralizing  encounter 
between  charge  carriers  of  opposite  sign  moving  about  in  a  gas  [Ref.  20]. 
The  first  studies  of  this  effect  were  made  by  Thompson  and  Rutherford  in 
1896.   They  stated  the  principles  of  recombination  to  explain  the  reduc- 
tion of  conductivity  in  a  partially  ionized  gas. 

The  recombination  coefficient  is  a  measure  of  the  number  of  recom- 
bining  events  per  unit  time  and  volume.  It  may  be  related  to  a  re- 
combination cross  section  [Ref.  21]. 

When  an  electron  and  an  ion  interact  or  collide  at  relatively  low 
velocity,  they  have  a  finite  probability  of  recombining  into  a  neutral 
atom.  In  order  to  conserve  momentum,  a  third  body  must  be  present. 
This  body  may  be  another  particle  or  photon.  If  it  is  a  photon,  then  it 
is  called  radiative  or  two-body  recombination.  When  another  particle  is 
involved,  it  is  called  three-body  recombination.  This  represents  a  loss 
term  to  the  continuity  equation  [Ref.  22].  Two-body  recombination 
predominates  when  the  gas  pressure  is  below  a  few  Torr. 

Due  to  losses  of  electrons  through  diffusion,  attachment  and  etc.  , 
an  accurate  value  of  the  recombination  coefficient  is  difficult  to 
obtain  experimentally.  Biondi  [Refs.  23  and  38]  has  found  that  the  time 
and  pressure  dependencies  of  electron  density  in  certain  gases  follow  a 
two-body  recombination  law  over  a  wide  range  of  variables.  The  gas  must 
be  quite  pure,  however,  to  reproduce  his  results  as  the  formation  of 
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negative  ions  and  the  Penning  effect  can  completely  destroy  any  sense  in 
the  experimental  data. 

Complications  arise  when  the  effects  of  diffusion  are  included  in 
the  recombination  measurements.  Some  of  these  effects  have  been  analy- 
zed for  infinitely  long  cylinders  [Ref.  24].  One  important  conclusion 
is  that  linearity  of  the  plot  of  one  over  the  electron  number  density 
versus  time  does  not  imply  a  meaningful  value  for  the  recombination 
coefficient  derived  from  the  slope  of  the  plot.  It  turns  out  that 
linearity  may  result  even  when  diffusion  is  controlling  the  decay 
process.  Either  the  diffusion  loss  term  or  the  recombination  loss  term 
may  be  dominant  over  most  of  the  plasma  at  any  particular  time  (except 
at  late  times  in  the  afterglow  when  diffusion  must  dominate  the  electron 
loss)  [Ref.  25].  Therefore,  the  plasma  may  be  either  diffusion  or 
recombination  controlled  at  any  given  time.  If  one  of  the  two  mechan- 
isms dominate  for  a  sufficiently  long  period  of  time  during  the  decay, 
then  the  electron  distribution  will  follow  a  solution  which  would  occur 
if  the  other  mechanism  were  absent.  (See  discussion  in  Chapters  V  and 
VI  on  linear  and  non-linear  solutions  of  the  electron  decay  equation.) 

J.  J.  Thompson  found  experimentally  that  diffusion  to  the  walls 
alone  could  not  account  for  the  reduction  of  conductivity  of  an  ionized 
gas.  The  recombination  of  a  free  electron  with  an  ion  was  a  significant 
factor  in  the  neutralization  process. 

Recombination  between  an  electron  and  an  ion  can  take  place  only 

_o 

when  the  electron  is  in  the  vicinity  of  the  ion  for  about  10  seconds. 
Clearly,  only  slow  moving  electrons  will  be  able  to  recombine  with-  the 
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ion  as  it  needs  to  remain  in  the  vicinity  of  the  ion.   During  two  body 

recombination,  the  electron  loses  excess  energy  through  photon  emission. 

Otherwise,  three  body  recombination  where  the  electron  transfers  excess 

energy  to  another  electron  is  the  dominant  process.  The  collision  cross 

section  for  the  electron-electron  collision  depends  on  the  electron 

temperature  and  is  thus  an  unlikely  event  in  normal  gases.  When  looking 

at  collision  times,  if  (t  )  is  the  smallest,  then  the  neutral  particles 

en-  r 

determine  the  recombination  rate  of  the  gas.  This  is  determined  purely 
by  species  pressure.  If  (t  .)  is  smallest,  then  forces  other  than 
species  pressure  determine  recombination  rates.  For  a  neutral  inter- 
action, the  collision  must  occur  within  the  Bohr  radius.  However,  for 
charged  particles,  the  Debye  rad" us  will  be  the  determining  factor. 
Since  the  Debye  radius  is  much  greater  than  the  Bohr  radius,  the  charged 
particle  interaction  will  dominate. 

When  ionizing  a  neutral  gas,  the  density  of  electrons  and  ions 
should  be  the  same  for  single  ionization.  The  number  of  recombinations 
per  second  is  then: 
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where  alfa  is  the  constant  of  proportionality.  Alpha  is  in  fact  not  a 
constant  [Ref.  26].  Consider  the  case  where  the  ionizing  source  is 
turned  off.  Many  of  the  free  electrons  will  still  be  near  their  former 
atom.  The  probability  of  electrons  recombining  with  the  ions  is  large. 
Therefore,  the  recombination  coefficient  will  initially  be  large  and 
then  become  smaller  with  time. 

It  at  time  zero,  the  density  of  charge  carriers  is  n  ,  one  can  inte- 
grate the  above  equation  directly.  This  becomes: 


(5) 


which  has  the  solution: 


n=  rS^Tt  <« 


This  is  called  the  ion  loss  formula  [Ref.  27]. 

Recent  investigations  show  that  the  number  of  charge  carriers  pro- 
duced by  ionization  may  be  different  due  to  attachment,  multiple  ioni- 
zation and  etc.  Also  it  has  been  found  that  two  different  ions  can 
exist  at  the  same  time  and  may  have  different  recombination  rates 
[Ref.  28]. 
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V.   LINEARIZED  BREAKDOWN  EQUATION  AND  SOLUTION 

A  model  of  the  streamer-breakdown  event  is  necessary  to  adequately 
describe  the  physical  events  which  occur.  The  streamer  is  thought  to  be 
a  channel  or  tube  of  high  charge  concentration  which  leads  to  high 
conductivity.  This  channel  of  ions  supports  the  flow  of  electrons  from 
the  cathode  to  the  anode  during  the  arc  transition.  The  effects  which 
reduce  the  ion  concentration  and  the  electron  density  have  already  been 
described  as  being  two  body  recombination,  diffusion,  convection  and 
electron  attachment.  Many  more  factors  need  to  be  included  in  any 
complete  description  of  the  events  that  take  place  within  the  streamer 
tube  (such  as  three-body  recombination). 

Radiative  recombination  (two-body  recombination)  requires  a  photon 
to  carry  away  excess  energy  from  an  electron-ion  recombination.  In 
three-body  recombination,  a  third  body  removes  excess  energy.  This 
third  body  may  either  be  an  electron  or  some  heavy  particle.  An  elec- 
tron may  also  recombine  with  a  molecular  ion.  If  the  excess  energy  goes 
into  tne  dissociation  of  the  molecule,  this  is  called  dissociative 
recombination.  Another  way  of  removing  excess  energy  is  to  form  a 
neutral  atom  in  which  two  electrons  are  simultaneously  excited.  This  is 
known  as  dielectronic  recombination. 

The  total  time  rate  of  change  of  the  charged  particle  density  may 
thus  be  written: 
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where  the  first  term  on  the  right  is  the  diffusion  term,  followed  by  the 
convective  term,  the  two  body  recombination  and  three  body  recombination 
terms.  (Da)>  (a?)  and  ^3)  are  the  respective  coefficients.  U  is  the 
flow  velocity  of  the  gas.  The  geometry  is  shown  in  Figure  2. 

First,  an  order  of  magnitude  test  is  made  to  determine  the  relative 
importance  of  these  terms.  To  do  this,  characteristic  times  are  defined 
and  their  magnitudes  compared.  Nominal  values  were  u:;ed  from  Table  1 
for  the  above  coefficients  for  nitrogen  at  atmospheric  pressures. 

Actual  values  of  the  above  coefficients  are  dependent  on  the  ionic 
species  in  the  gas.  It  is  still  possible  to  obtain  general  order  of 
magnitudes  for  the  coefficients,  however.  The  values  shown  were  taken 
from  Mitchner  and  Kruger  [Ref.  30].  Because  (Von>>t  ^  tne  convectlve 
term  is  neglected. 

These  characteristic  times  show  that  the  effect  of  recombination  is 
considerably  faster  than  all  but  turbulent  flow  on  the  alteration  of 
charged  particle  density.  For  a  first  order  approximation,  then,  one 
could  assume  negligible  convective  and  diffusion  effects  if  the  fluid 
flow  were  laminar.  This  has  already  been  done  without  stating  so  in  the 
chapter  on  recombination.  There  the  solution  was  found  to  be: 
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Figure  2.   Streamer  Geometry 
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Table  1.   Coefficient  Values  for  Nitrogen 
[Ref.  29] 


No 

— 

1020 

,r3 

aZ 

- 

lO"13 

3/sec 

a3 

- 

10"35 

3/sec 

ro 

= 

10"4 

m 

Da 

= 

lO'4 

m2/sec 

for 

laminar  flow 

Da 

3 

10"1 

2/sec 

rn 

for 

turbulent  flow 

L 

= 

10"1 

m 

u 

= 

102 

in/sec 

T 

conv 

= 

L/U 

=  lO'3 

sec 

Tdiff 

= 

R  2/D 
o   a 

=  10"4 

sec 

for  laminar  flow 

Tdiff 

= 

R  2/D 

o   a 

=  10"7 

sec 

for  turbulent  flow 

i 
recom 

- 

1/ttN 

0 

=  10"7 

sec 
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niti  = 


J  +  ^lOot  (8) 


with  the  boundary  conitions  of: 


fl(o)  =  tf„    A*4    fl(«D  =  O  (9) 


If  the  flow  is  fully  turbulent,  the  diffusion  term  may  not  be 
neglected  in  the  calculations.  Therefore,  let  us  assume  that  the  ga-s 
flow  is  in  some  intermediate  state;  neither  laminar  nor  turbulent. 
Turbulent  fluid  flow  may  increase  the  diffusion  coefficient  by  three  or 
more  orders  of  magnitude.  The  characteristic  diffusion  time  will  then 
be  on  the  same  order  of  magnitude  as  the  recombination  time. 

The  equation  which  models  the  situation  for  two-body  recombination 
is: 


If-   ^£«?)— ' 


r 


The  three-body  analysis  is  given  in  Appendix  A.   To  achieve  a  linear 
approximation  of  this  equation  we  must  insure  the  flow  is  only  mildly 
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turbulent.  With  this  in  mind,  diffusion  may  be  introduced  as  a  "pertur- 
bation" of  the  recombination  equation.  Equation  (10)  can  then  be 
linearized  because  if  the  perturbation  is  small,  then  the  non-linear 
terms  can  be  neglected.  First  the  charge  density  is  written  in  two 
parts:  a  base  and  a  perturbation  part.  The  base  charge  density  ex- 
presses the  state  of  the  plasma  channel  in  the  absence  of  the  diffusion. 
The  linear  theory  is  valid  as  long  as  the  diffusion  term  is  much  smaller 
than  the  recombination  term.  This  is  true  for  non-turbulent  situations. 
To  induce  the  fully  turbulent  case,  a  numerical  solution  of  the  govern- 
ing equation  must  be  done.  This  will  follow  in  Chapter  VI. 
The  charge  carrier  density  is  now  written  in  the  form: 


rUr.-t)  -  V^CtW  i^Cf.t)  lA>Uere_   £  «  i_ 


(11) 


In  order  to  get  an  analytic  solution  and  make  the  derivation  easier 
to  ^ollow,  we  shall  switch  from  the  axisymmetric  or  polar  case  to  a  two 
dimensional  solution.  The  governing  equation  is  now: 


at  oxfc  (12) 
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where  now  we  let: 


ni*.*!*  Xu\  +  cn,(%,iz) 


(13) 


This  is  substituted  into  the  governing  equation: 


-    -ot^n%.at^jo,  +  2VV] 


^  fl.AT. 


Z   J        d5) 


Now  we  may  neglect  the  higher  order  terms  and  also  cancel  the  base 
part: 

2t  Ox3-  (16) 


JfolKJot-  (17) 
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The  linearized  solution  to  equation  (10)  is: 


«1(% 


,or  *<>■*  ♦  £* 


^ft&MftftS 


i+*«.t    i(IWwrf  L      viTSi  j f  cr*  \  wSvJ  (is) 

•for   X±X* 


&?  fc<o»)  ♦«*(»)] 


(18a) 


X>Xo 


The  error  function  correctly  models  the  step  decrease  in  the  elec- 
tron density  at  the  boundary  at  time  zero.  The  utility  of  the  error 
function  is  also  demonstrated  at  times  greater  than  zero.  The  function 
allows  the  effusion  of  electrons  across  the  boundary  of  the  ionized 
channel.  This  process  models  the  effect  of  diffusion  while  still 
meeting  the  requisite  boundary  conditions.  A  plot  of  the  error  function 
in  non-dimensional  parameters  is  shown  in  Figure  3.  The  three-body 
counterpart  of  equation  (18)  is  given  in  Appendix  A. 

The  diffusion  coefficient  appears  only  in  the  argument  of  the  error 
function.  The  effect  of  larger  diffusion  coefficients  is  to  increase 
the  rate  at  which  the  electrons  leave  or  diffuse  from  the  original 
plasma  channel.  If  the  diffusion  coefficient  were  zero,  then  the  elec- 
trons would  never  cross  the  initial  boundary  of  the  plasma  channel. 
That  is,  the  form  of  the  channel  would  be  independent  of  time.  However, 
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Figure  3.   Plot  of  Error  Function  (ERF) 
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the  electron  density  within  the  channel  would  decrease  in  time  due  to 
the  recombination  term.  These  numerical  results  are  in  keeping  with  the 
physical  picture  of  the  plasma  channel. 

It  must  be  remembered  that  all  discussion  pertaining  to  the  dif- 
fusion of  electrons  really  means  the  ambipolar  diffusion  of  the  elec- 
trons and  ions  together.   This  was  described  in  Chapter  III. 
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VI.   THE  NON-LINEAR  EQUATION  AND  SOLUTION 

Non-linear  equations  are  defined  as  those  which  contain  products  of 
the  dependent  variables  or  transcendental  functions.  One  can  solve 
these  equations  either  for  the  "roots"  or  for  solutions  which  satisfy  a 
set  of  non-linear  equations.  The  solution  to  equation  (7)  which  models 
the  de-ionization  of  the  plasma  channel  seeks  the  roots. 

There  appears  to  be  no  solution  in  terms  of  known  functions  to  the 
general  form  of  equation  (7).  Sometimes  it  is  not  clear  whether  a  given 
equation  has  a  solution  at  all.  The  solution  of  the  non-linear  form  is 
of  an  interative  nature.  To  ensure  convergence,  it  was  necessary  that  a 
reasonable  approximation  to  the  solution  in  question  be  known  prior  to 
iterating  [Ref.  31].  That  was  the  purpose  of  the  linearized  solution 
presented  in  the  precaedirg  chapter.  If  the  solutions  at  all  match, 
then  there  may  be  some  degree  of  confidence  that  the  solutions  are 
correctly  describing  the  situation.  If  there  is  a  great  disparity 
between  the  two  solutions,  then  the  reason  for  the  discrepancy  must  be 
resolved. 

Because  the  geometry  of  a  streamer  is  cylindrical  in  nature,  the 
cylindrical  coordinate  system  is  used  in  the  numerical  analysis.  The 
system  is  then  described  by  two  dimensions;  radial  position  and  time. 
It  is  assumed  that  there  exists  angular  symmetry  and  that  the  length  of 
the  ionized  channel  is  much  greater  than  it's  diameter.  This  permits 
the  exclusion  of  electrode  effects  at  either  end. 
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The  specification  of  the  boundary  conditions  was  not  without  some 
ambiguity.  Boundary  conditions  are  specified  along  the  centreline  of 
the  channel  and  at  infinity  with  a  finite-difference  computational 
method. 

In  the  computational  model,  the  two-dimensional  axes  represent  tima 
in  the  vertical  and  radial  position  along  the  horizontal.  At  time  zero, 
the  electron  density  in  the  radial  direction  is  presumed  to  be  constant 
out  to  some  point  r  .   After  this  point,  the  electron  density  Is  zero. 

Along  the  centreline  of  the  plasma  channel,  the  slope  of  the  elec- 
tron density  versus  position  must  be  zero.  This  is  used  to  define  the 
starting  value  for  the  electron  density  along  the  centreline  for  any 
time  greater  than  zero. 

Knowing  these  boundary  conditions,  the  geometry  of  the  situation 

lends  itself  well  to  a  finite  difference  computational  molecule.   The 

grid  size  in  the  radial  direction  is  relatively  insensitive  in  the 

solution  of  the  equation.   A  convenient  step  size  of  r  /10  is  used.   The 

same  is  not  true  for  the  step  size  in  the  time  dimension.   The  primary 

reason  for  this  leads  back  to  the  characteristic  times  encountered  in 

the  linearized  solution.   The  timewise  step  size  needs  to  be  smaller 

than  the  smallest  characteristic  time  of  the  system.   In  this  case  it 

turned  out  to  be  the  recombination  time.   The  timewise  step  size  was 

thus  taken  to  be  t   /10. 

rec 

The  technique  used  along  with  the  computational  molecule  is  the 
explicit  method  of  temporal  development  [Ref.  32].  The  principal  defect 
in  this  form  of  solution  is  the  need  for  very  small  time  increments  to 
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ensure  stability.  This  leads  to  lengthy  computation.  A  reasonable 
modification  to  this  solution  to  extend  it  to  greater  times  without  the 
need  for  many  iterations  would  be  to  go  to  an  implicit  recurrence  formu- 
lation. This  is  where  two  or  more  unknown  values  of  the  equation  at 
t(i+l)  are  expressed  in  terms  of  known  values.  This  can  be  done  by  the 
solution  of  sets  of  simultaneous  equations,  which  complicate  the  solu- 
tion considerably. 

What  was  originally  deemed  a  "stability  border"  was  found  in  both 
the  linear  and  non-linear  solutions.  This  border  governed  the  size  of 
the  time  steps.  If  the  delta  time  was  too  large,  solutions  which  con- 
tained physically  meaningless  densities  (negative  values)  were  the 
result.  This  border  was  found  to  be  roughly  one-fifth  to  one-tenth  the 
shortest  characteristic  time  in  the  system.  This  empirical  finding  lead 
to  ';he  development  of  the  one-tenth  of  the  recombination  time  as  the 
time  step  interval.  While  researching  the  literature,  this  small  time 
step  was  shown  to  be  a  necessary  consequence  of  the  numerical  method 
employed  as  previously  stated. 

If  one  considers  a  cartesian  coordinate  system  where  the  'X'  axis 
represents  radial  position  and  the  ' Y'  axis  is  time,  the  values  for  the 
electron  density  can  be  established  using  a  "marching"  process  [Ref.  33]. 
The  boundary  conditions  specify  the  values  of  "n"  (electron  density) 
along  the  radial  position  at  time  zero.  The  requirement  of  zero  slope 
of  the  density  curve  at  the  centre  of  the  plasma  'tube'  for  any  time  is 
a  second  requirement.  Thus,  n(r,0)  is  specified.  The  density  at  other 
points  is  then  determined  by  the  density  at  preceeding  points. 
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With  the  marching  scheme  in  mind,  the  numerical  evaluation  of  the 
points  on  the  graph  are  done  in  the  following  manner: 

i  =  radial  position  marker 

j  =  time  increment  marker 


DtJiJ    At  1  ^.a-Aij 


(19) 


+&(r$*<fc*hj  i******]****  K<  ->w  \  (20) 


There  is  a  restriction  that  is  easily  met  in  this  solution  where: 


Ar*  >  iDa  At 


[Ref.  31]    (21) 


Because  of  the  small  time  steps  used,  this  restriction  has  no  effect  on 
the  calculations  unless  the  user  reduces  the  radial  step  size  to  very 
small  steps.   In  the  physical  picture,  the  initial  radius  of  the  dis- 
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charge  is  on  the  order  of  a  few  Debye  lengths.  However,  to  reduce  the 
step  size  much  smaller  than  one-tenth  to  one-hundredth  of  the  discharge 
radius,  the  time  interval  will  have  to  be  adjusted  accordingly. 

These  two  parts  of  the  non-linear  equation  are  the  only  ones  that 
need  modification  to  perform  the  numerical  solution.  The  recombination 
terms  need  no  adjustments  as  there  is  no  time  or  positional  derivatives 
associated  with  them.  They  depend  principally  on  the  charge  density  of 
the  system. 

The  governing  equation  which  is  originally  written  as: 


Dt   "r  0<-lr  Or*;  -  -oLn 


(22) 


may  now  be  written: 


Lt. 


fc**M^N"^*H"&l^-M"-**<»> 


This  is  the  form  used  in  the  numerical  solution.  The  program  de- 
fines the  initial  conditions  and  boundary  conditions  before  iteration 
begins.  It  then  progresses  along  the  radial  position  axis  at  a  fixed 
time  increment  based  upon  the  shortest  characteristic  time  in  the  system 
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and  evaluates  the  charge  density  at  every  point  and  stores  it  in  an 
array.  When  the  program  has  determined  the  charge  density  at  all  points 
in  the  array,  it  then  proceeds  to  plot  the  densities  in  one  of  several 
outputs  specified. 

All  constants  in  the  program  may  be  changed.  The  only  care  that 
need  be  observed  is  when  adjusting  any  parameter  which  changes  the 
characteristic  recombination  or  diffusion  times.  The  time-wise  iter- 
ation steps  are  based  on  the  recombination  time  as  is  the  diffusion 
time.  Thus,  in  general,  no  modification  to  the  program  is  required  when 
changing  initial  charge  c'ensity,  n  ,  or  the  recombination  coefficient,  a. 

If  one  wishes  to  make  the  diffusion  time  shorter  than  the  recombina- 
tion time,  the  time-wise  step  will  need  to  be  modified  such  that  it  is 
approximately  one-tenth  of  the  shortest  characteristic  time  of  the 
system. 
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VII.   EFFECT  OF  THE  DIFFUSION  COEFICIENT 

Experiments  have  shown  that  the  diffusion  coefficient  changes  its 
value  between  laminar  and  turbulent  flow.  The  difference  is  an  incraase 
of  approximately  three  orders  of  magnitude.  This  can  reduce  the  char- 
acteristic diffusion  times  under  turbulent  flow  to  the  same  order  as  the 
recombination  time.  Turbulent  diffusion,  therefore,  could  not  be 
modeled  by  the  linearized  solution.  As  shown  in  Appendix  B,  the  dif- 
fusion times  in  the  linearized  cases  were  allowed  to  vary  between  ten- 
thousand  to  ten  times  the  recombination  times.  To  summarize  those 
results,  one  can  say  that  in  all  cases,  the  charge  density  diffused 
outward  from  the  plasma  channel  at  a  greater  rate  when  the  diffusion 
coefficient  increased  from  laminar  to  turbulent. 

Turbulence  has  a  strong  influence  on  the  characteristics  of  elec- 
trical discharges.  This  is  mainly  due  to  two  effects.  First,  the 
ambipolar  diffusion,  which  also  determines  the  rate  of  loss  of  charged 
particles  to  the  walls,  is  greatly  enhanced.  In  diffusion  dominated, 
self-sustaining  plasma  columns,  the  balance  is  thus  disturbed  and  this 
causes  a  substantial  increase  in  the  ionization  rate,  local  electric 
fields  and  electron  temperature. 

Boundary  layers  seperate  and  mix  the  flow,  thus  increasing  the 
transport  coefficients  of  the  plasma,  in  particular  its  heat  conduc- 
tivity and  ambipolar  diffusion.  The  increase  in  ambipolar  diffusion 
causes  the  electron  density  profile  to  be  more  uniform  over  the  cross 
section  [Ref.  34]. 
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Turbulent  gas  flow  means  that  at  every  point  in  the  gas  the  velocity 
fluctuates.  That  is,  the  velocity  changes  its  magnitude  and  direction. 
When  the  velocity  is  steady,  there  are  three  fluctuating  velocity  com- 
ponents whose  time  average  vanishes  over  sufficiently  long  periods  of 
time. 

Turbulence  involves  statistical  fluctuations  and  thus,  no  two  tur- 
bulent streams  are  alike.  There  are  two  average  properties  which  are 
used  to  describe  turbulence.  These  are  the  turbulent  intensity  and  the 
length  or  scale  of  turbulence.  When  the  scale  of  turbulence  is  small 
compared  to  the  dimensions  of  the  system  involved  (here  it  is  the 
streamer  radius),  the  turbulent  intensity  alone  suffices  to  characterize 
the  flow  [Ref.  35].  This  is  useful  in  hydrodynamics  and  aerodynamics. 
rhis  simplification  is  not  valid  in  the  glow-discharge  situation  because 
both  glow-discharge  and  turbulence  are  microscopic  properties.  This 
tends  to  make  precise  analysis  difficult  at  best. 

As  the  turbulence  of  the  free  stream  increases,  an  increase  in  the 
rate  of  heat  transfer  results  due  to  the  flow  of  mass  from  the  system 
under  investigation.  The  effect  is  remarkably  great  as  a  turbulent 
intensity  of  approximately  2.5%  produces  an  increase  in  the  local  heat 
flux  of  about  80%  [Ref.  36]. 

The  question  which  may  be  asked  here  is  whether  the  charge  would 
recombine  at  a  sufficiently  rapid  rate  so  as  to  reduce  the  charge 
density  before  it  has  a  chance  to  diffuse  out  of  the  plasma  channel. 
This  is  a  result  of  the  linear  solution  assumptions  and  the  evaluation 
of  the  characteristic  time  of  diffusion  and  recombination.  If  the 
recombination  times  are  in  fact  several  orders  of  magnitude  faster  than 
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diffusion,  then  the  ionization  density  within  the  plasma  channel  will 
decrease  rapidly  before  the  charge  carriers  can  diffuse.  This  will  mean 
that  diffusion  is  really  not  the  dominant  charge  loss  mechanism.  This 
turns  out  to  not  be  the  case  at  all.  In  both  the  linear  and  non-linear 
solutions,  the  charge  carriers  diffuse  out  of  the  plasma  channel  long 
before  the  density  drops  appreciably  due  to  recombination. 

It  has  been  previously  noted  in  both  Chapters  II  and  III  that  dif- 
fusion is  the  dominant  loss  mechanism  in  the  problem  posed.  This  has 
been  verified  experimentally  for  some  time  [Ref.  37].  However,  one  must 
focus  on  the  assumption  that  at  time  zero,  the  charge  density  is  a  step 
function;  that  is,  the  charge  density  is  zero  outside  the  plasma  channel. 
Originally,  a  very  steep  density  gradient  (on  the  order  of  infinity) 
exists.  This  acts  as  a  very  strong  drivirg  mechanism  for  diffusion. 
The  effects  of  diffusion  at  very  short  time  intervals,  therefore, 
dominate  over  recombination. 

The  results  in  Figures  4  through  9  show  that  diffusion  is  in  fact 
the  dominant  charge  loss  mechanism.  The  figures  clearly  demonstrate 
this  for  all  but  a  highly  laminar  flow  situation.  The  figures  all  plot 
the  log(n)  versus  radial  position.  The  temporal  development  of  the 
charge  density  profile  is  then  shown  as  a  family  of  curves.  The  charge 
loss  due  to  recombination  is  somewhat  lost  due  to  the  logrithmic  func- 
tion. This  was  necessary  to  show  clearly  both  the  reduction  in  charge 
density  and  its  diffusion  out  of  the  plasma  channel. 
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Figure  4.   Results  of  Non-Linear  Program  for  p  =  10 
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Figure  5.   Results  of  Non-Linear  Program  for  0  =  5,000 
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Figure  6.   Results  of  Non-Linear  Program  for  p  =  1,000 
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Figure  7.   Results  of  Non-Linear  Program  for  0  =  500 
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Figure  8.   Results  of  Non-Linear  Program  for  0  =  100 
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Figure  9.      Results  of  Non-Linear  Program  for  p  =  50 
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Table  2  lists  the  plotted  character  and  the  associated  time  it 

represents  on  the  preceeding  figures.   0  is  the  ratio  of  xn/TR.   It  is 

3 
characteristic  of  the  degree  of  turbulence  in  the  flow.   When  p  >  10  , 

the  flow  is  considered  to  be  laminar.   If  0  <  50,  then  the  flow  is 

considered  to  be  fully  turbulent. 


Table  2.  Plotted  Character  Definitions 
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The  graphs  on  the  preceeding  pages  are  the  results  of  the  full 
non-linear  equations  as  analyzed  in  the  computer.  The  technique  has 
been  outlined  in  Chapter  VI  and  the  program  used  is  shown  in  Appendix  B. 
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All  of  the  plots  show  the  temporal  development  of  the  plasma  cross- 
section.  This  is  achieved  by  graphing  the  charge  density  versus  radial 
position.  The  family  of  curves  on  each  plot  is  then  the  charge  density 
profile  at  times  greater  than  zero.  Figure  4  is  for  a  laminar  diffusion 
coefficient,  Figure  9  is  for  a  turbulent  diffusion  coefficient,  and 
Figure  5  through  Figure  8  are  for  arbritary  intermediate  levels  of 
diffusion. 

Figure  4  shows  the  temporal  development  of  the  charge  density  for  a 
laminar  gas  flow.  Note  the  very  steep  density  gradient  at  time  zero. 
This  gradient  acts  to  drive  the  diffusion  process  rapidly  at  first. 
However,  as  the  gradient  is  reduced  in  time,  the  driving  force  is  also 
reduced.  Because  of  this,  the  charge  density  profile  remains  very  much 
like  that  of  the  original  plasma  channel. 

This  relatively  stable  region  of  high  charge  density  translates  to  a 
stable  region  of  low  resistance  to  the  flow  of  electrons.  Laoking  at 
the  situation  physically,  future  breakdown  streamers  will  be  enticed  to 
follow  this  highly  conductive  path  and  either  complete  the  breakdown 
process  if  it  has  not  already  occurred,  or  they  will  sustain  the  break- 
down. This  stable  channel  will  permit  breakdown  to  occur  at  lower 
current  levels  once  the  channel  has  been  established. 

Figure  9  is  at  the  other  end  of  the  diffusion  spectrum,  that  of  a 
turbulent  gas  flow.  Here,  the  diffusion  coefficient  is  several  orders 
of  magnitude  larger  than  in  the  laminar  flow  regime.  The  associated 
graph  clearly  demonstrates  the  rapid  increase  in  the  plasma  channel's 
cross-section  with  time.   This  growth  acts  to  reduce  the  overall  charge 
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density  within  the  plasma  channel.  Future  streamers  will  see  a  higher 
path  resistance  than  in  the  laminar  flow  case.  If  the  path  were  created 
by  a  streamer  which  did  not  cause  breakdown,  future  streamers  will  be 
less  likely  to  follow  the  path  and  complete  the  breakdown  process.  If 
the  plasma  channel  is  a  result  of  the  breakdown  process,  future  stream- 
ers will  need  to  be  more  energetic  to  sustain  the  breakdown  if  they 
follow  the  original  path. 

As  a  result  of  the  ionization  process  and  the  conclusions  derived 
from  the  analysis  of  the  'governing  equation1,  it  is  hypothesized  that 
turbulent  flow  of  the  gaseous  medium  in  a  glow  discharge  process  will 
delay  the  glow  discharge  to  arc  transition  when  compared  to  a  laminar 
flow  of  gas. 
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VIII.   CONCLUSIONS 

This  thesis  presents  the  background  concerning  the  process  of  de- 
ionization  by  two-body  recombination  and  ambipolar  diffusion  in  a  gas. 
From  this  background,  a  description  of  the  system  is  developed  and 
analyzed.  The  results  from  this  analysis  demonstrate  the  effects  of 
diffusion  when  the  flow  is  either  laminar  or  turbulent. 

These  results  show  that  with  the  increase  in  turbulence  and  thus 
with  the  increase  in  the  diffusion  coefficient  the  charge  density 
profile  spreads  out  rapidly.  This  radial  expansion  reduces  the  charge 
density  within  the  streamer  discharge  'tube".  The  results  are  somewhat 
inconclusive,  however,  when  considering  the  physical  process  which  takes 
place  afterwards.  To  explain  more  clearly  the  question  at  hand,  an 
examination  of  the  post-streamer  events  in  two  limiting  cases  must  first 
be  examined. 

First,  assume  the  pure  laminar  flow  situation  whether  there  is 
negligible  diffusion.  Here,  the  streamer  will  retain  its  initial  charge 
density  profile  as  the  charge  carriers  are  lost  by  recombination.  The 
streamer  will  not  grow  in  diameter,  but  the  charge  density  will  decay 
more  slowly  than  when  diffusion  is  added.  This  streamer  "tube"  will 
maintain  its  high  conductivity  as  the  conductivity  of  a  partially  ion- 
ized plasma  is  proportional  to  the  charge  density. 

The  second  case  deals  with  a  highly  turbulent  flow  where  there  exist 
very  large  diffusion  effects.   The  streamer  in  this  situation  will  have 
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a  rapidly  expanding  charge  density  profile.  The  charge  carriers  within 
the  streamer  will  still  be  lost  through  recombination  and  the  charge 
density  will  be  further  reduced  through  diffusion. 

In  the  limit  of  numerous  streamer  events,  the  second  case  will 
establish  a  more  uniform  distribution  of  charge  density  throughout  the 
gas.  However,  the  first  case  of  no  diffusion  will  have  numerous 
streamer  "tubes"  propagating  through  the  gas,  retaining  their  original 
charge  density  profile  and  decaying  more  slowly  due  to  recombination. 

Now,  the  question  of  follow-on  streamers  may  be  examined.  Will  the 
later  streamers  seek  out  the  narrow,  highly  conductive  channels  in  which 
to  propagate  further  or  will  they  propagate  independently  until  the 
charge  density  within  the  gas  as  a  whole  is  large  enough  to  sustain  a 
breakdown  transition?  Intuitively,  either  case  may  be  argued  success- 
fully to  some  extent.  Further  research,  more  than  likely  experimental, 
is  necessary. 

In  addition  to  the  follow-on  streamer  problem,  another  interesting 
process  might  be  worthy  of  consideration.  Given  the  case  of  near 
turbulent  flow,  where  diffusion  acts  to  change  the  charge  density  pro- 
file, what  takes  place  at  relatively  long  time  periods  within  the 
expanding  streamer  "tube"?  Remember  that  the  recombination  rate  is 
proportional  to  the  square  of  the  charge  density  (the  number  of  ions 
times  the  number  of  free  electrons).  If  the  recombination  coefficient 
is  constant  throughout  the  streamer  "tube",  areas  of  high  density  will 
decay  faster  than  areas  of  lower  density  (as  in  the  diffusion  expansion 
region).   This  process  could  indicate  the  existence  of  a  perturbation  in 
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the  charge  density  profile  where  the  charge  density  is  either  less 
inside  this  location  and  radially  outward  or  greater.  The  result  would 
be  some  type  of  higher  or  lower  density  region  (saddle  point)  at  a  given 

radial  position  within  the  charge  density  profile. 

-5      -4 
If  this  process  is  correct  for  times  on  the  order  of  10  '  to  10 

seconds,  then  there  may  be  an  optimum  level  of  diffusion  which  delays 

the  process  of  streamer  to  arc  breakdown  or  enhances  the  transition. 

These  questions  remain  unresolved  at  this  point.   Further  research 

into  may  lead  tc  an  optimum  level  of  turbulent  mixing  within  the  gaseous 

medium  to  maximize  the  ionization  level  within  the  gas  and  yet  not  allow 

the  glow  discharge  to  arc  transition.   This  would  maximize  either  the 

power  output  or  the  safety  characteristics  of  electrically  pumped  lasers. 
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APPENDIX  A 
THREE-BODY  LINEAR  SOLUTION 

The  following  is  the  mathematical  derivation  of  the  three-body 
governing  equation  linearization  process.  Results  are  given  in  terms  of 
the  Gaussian  Error  Function  (erf)  for  the  analytical  solution  of  the 
two-dimensional  case. 

The  non- linear  equation  which  describes  the  continuity  equation  for 
charge  density  in  the  cartesian  coordinate  system  is: 


3^-^*—  *-±J\  t*     +^0-bod1  (Al) 


(A2) 


The  assumptions  required  are: 

(1)  Diffusion  is  a  perturbation  on  the  recombination  controlled 
process 

(2)  xR  «  xD 
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Three-body  linearized  equation  derivation: 


D-b  D>*  (A3) 


*2»     -±*J 
D-b 


yields : 


=    -J-[*V3tiiln',>  3Sj^n  nib? J 
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(A4) 


"'     L^a]"  (A5) 


let:       AUtt^a      Xli)f  f  «,M      u>U*r<.      €  4.<  1  (A6) 


(A7) 


(A8) 


because  of  equation   (A4): 


(A9) 


2*i-U.  SIS  =   -3dn.n/- 

Ot  O**-  (A10) 


When  the  appropriate  boundary  conditions  are  included  the  specific 
solution  is: 


AHA*™*    .♦  12 


(All) 
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APPENDIX  B 
NGN-LINEAR  PROGRAM 

The  program  listing  shows  several  features  which  need  amplification. 
The  first  pertains  to  the  centreline  condition  and  how  it  effects  the 
computational  molecule.  One  bourdary  condition  requires  the  slope  of 
the  density  profile  be  zero  at  the  centreline.  This  is  used  in  the 
program  to  evaluate  the  charge  density  at  the  station  next  to  the  cen- 
treline. This  value  is  the  same  as  the  centreline  value.  Due  to  this 
technique,  the  computational  molecule  does  not  need  to  "reach  across" 
the  centreline  to  find  the  (0-1, j)  magnitude.  The  program  uses  the 
value  at  (0,j)  and  at  (2,j)  to  evaluate  the  density  at  (l,j+l). 

Figures  4  through  9  in  Chapter  VII  are  plotted  to  five  times  the 
initial  streamer  radius.  However,  the  program  calculates  the  charge 
density  profile  out  to  nineteen  times  the  initial  radius  for  each  com- 
putation. The  region  not  shown  in  the  figures  may  be  plotted  by  suit- 
able changes  to  the  plot  commands  in  the  program. 

In  the  program,  the  diffusion  coefficient  (D  )  is  based  on  the 

a 

recombination  time  and  the  factor  p.  As  p  changes  from  laminar  to 
turbulent  conditions,  D  is  changed  accordingly.  This  technique  was 
utilized  to  quickly  adjust  the  flow  conditions  in  the  program. 

As  a  test  of  the  stability  condition  shown  in  equation  (21),  the 
following  worst  case  situation  is  presented.  The  largest  value  of  D  is 

a 

governed  by  the  choice  of  p,  as  the  other  parameters  are  constant  during 
iteration. 


69 


t>a  = 


Jivi.r,* 


I3 


Substituting  the  smallest  value;  p  =  50  yields  a  D  of: 

a 


J-flo^  ^    r3 


J-    =     S**/o 


brx  >£t>*  fct 


2.5-%/d"'    >   Z*to" 


(Bl) 


(32) 


This     is    the    largest    value    of    D      possible    in    this    program.       If    this 

a 

maximum  D  is  substituted  into  equation  (21)  we  find: 

3 


(B3) 


(B4) 


This  confirms  the  assumed  stability  of  the  solution  technique  for 
values  of  (3  >  50.  One  could  reverse  the  process  and  find  the  smallest 
value  for  p  which  satisfies  the  stability  criteria  if  necessary. 
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FILE:  NONLIN   FORTRAN   Al  NAVAL  POSTGRADUATE  SCHOOL 

//WALLACE  JOB  (2246,0096),'WALLACE,R,JfCLASS=A 

//EXEC  FRTXCLGP 

//FORT. SYS IN   DD   * 

C 

C 

C       THE  NON-LINEAR  SOLUTION  FOR  NEAR  TURBULENT  FLOW 

C 

C 

REAL  R0,N0,DA,ALFA1,ALFA2,DELTAT,DELTAR,N(199,199),A,B,C,INT,TIME 

REAL  VAL(199),RR(199) 

INTEGER  I, J, K,L, COUNT, KEY 
C 
C 

C        INITIALIZE  THE  CONSTANTS 
C 
C 

R0-50.E-6 

N0-1.E20 

ALFA1-1.E-12 

ALFA2-1.E-34 

DELTAT-1. / (alf al*N0*10. ) 

DKLTAR-RO/10.- 
C 
C 

C       CALCULATE  THE  DIFFUSION  COEFFICIENT  BASED  ON  RECOMBINATION 
C  TIME  AND  BETA 

C 
C 

DA-(ALFAl*N0*R0*R0)/50. 
C 

c 

C       IN  THE  ABOVE  EQUATION,  DIVIDE  BY  5000  FOR  LAMINAR  FLOW 

C  AND  3Y  50  FOR  TURBULENT  FLOW 

C 

£       SET  THE  MATRIX  TO  VIRTUAL  ZERO 

C 

C 

DO  100  J-1,199 
DO  200  1=1,199 
N(I,J)=1. 
200    CONTINUE 
100  CONTINUE 

DO  10  1-1,11 
N(I,1)=N0 
10  CONTINUE 

DO  300  J-1,199 

TIME-DELTAT*J 
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DO  400  1=2,191 
R=(I-1)*R0/10. 

A=(N((I+1),J)-(2.*N(I,J))+N((I-1),J))/(DELTAR*DELTAR) 
B=(N((I-i-l),J)-N((I-l),J))/(R*2.*DELTAR) 
C=(A*(A+B)-(ALFA1*N(I,J)*N(I,J))-(ALFA2*N(I,J)**3.) 
INT=R/R0 
IF  (I.NE.2)  GO  TO  111 

N((I-1),(J+1))=N(I,(J+1)) 
111    CONTINUE 
400  CONTINUE 
300  CONTINUE 
C 

C         TRANSFER  THE  RESULTS  FOR  THE  PLOT  ROUTINE 
C 

DO  8478  J-1,191 

DO  8477  1=1,50 

VAL(I)=(LOG(N(I,J)+1.))/2.303 

8477  CONTIN&E(I)=(I-1)/10- 
C 

C         VERSATEK  PLOT  ROUTINE 
C 

8478  CONTINUE 

CALL  PLOTG(RR,VAL,50,L,1,L, 'R/RO',4 , "LOG(N) ' ,6,0,0,0,0,6,6) 
CALL  PLOT (0.0, 0.0, 999) 
999  FORMATC  ' , 29X,F6. 2, 10X.E10. 3) 
888  FORMATC  ') 
555  FORMATC  ' ,E9.2, 5X,E9 .2) 
STOP 
END 
/* 

//GO.SYSIN  DD   * 
/* 
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